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Abstract — We are given an equal number of mobile robotic 
agents, and distinct target locations. Each agent has simple 
integrator dynamics, a limited communication range, and 
knowledge of the position of every target. We address the 
problem of designing a distributed algorithm that allows the 
group of agents to divide the targets among themselves and, 
simultaneously, leads each agent to reach its unique target. 
We do not require connectivity of the communication graph at 
any time. We introduce a novel assignment-based algorithm 
with the following features: initial assignments and robot 
motions follow a greedy rule, and distributed refinements of the 
assignment exploit an implicit circular ordering of the targets. 
We prove correctness of the algorithm, and give worst-case 
asymptotic bounds on the time to complete the assignment as 
the environment grows with the number of agents. We show that 
among a certain class of distributed algorithms, our algorithm 
is asymptotically optimal. The analysis utilizes results on the 
Euclidean traveling salesperson problem. 

I. Introduction 

CONSIDER a group of n mobile robotic agents and 
n target locations, all lying in R d , d > 1. Each 
agent has a limited communication range, and knows the 
location of some subset (possibly all) of the n targets through 
GPS coordinates or a map of the environment. The target 
assignment problem we consider is to design a distributed 
algorithm that allows the group of agents to divide the n 
targets among themselves and, simultaneously, that leads 
each agent to reach its unique target. Such a problem could 
arise in several applications. For example, one could think 
of the agents as UAV's on a surveillance mission, and the 
targets as the centers of their desired loitering patterns. Or in 
the context of formation control, the target positions could 
describe the desired formation for a group of robots. 

The first question is; how do we divide the targets among 
the agents in a centralized fashion? A reasonable strategy 
would be to minimize the sum of the distances traveled by 
each agent to arrive at its target. The problem of optimally 
dividing n persons among n objects, subject to a linear cost 
function, is a problem in combinatorial optimization [1]. It is 
referred to as the assignment problem, or the minimum weight 
perfect matching problem in bipartite graphs. The assign- 
ment problem can be written as an integer linear program. 
Unlike some integer linear programs, such as the Euclidean 
traveling salesperson problem (ETSP), optimal solutions for 
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the assignment problem can be computed in polynomial time. 
In 1955 Kuhn [2] developed the Hungarian method — the 
first polynomial solution for the assignment problem. Kuhn's 
method solves the problem in 0(n 3 ) computation time (see 
Section for a definition of the O notation). 

Another approach to the assignment problem is the auction 
algorithm [3], [4], [5], first proposed by Bertsekas. This 
method solves the problem in 0(n 3 ) computation time, but 
can be computed in a parallel fashion, with one processor 
for each person. Recently, Moore and Passino [6] modified 
the auction algorithm to assign mobile robots to spatially 
distributed tasks in the presence of communication delays. 
However, in order to exchange bids on a particular object 
(task), the auction algorithm, and thus the work in [6], 
requires that the communication graph between processors 
(robots) is complete. 

In this paper we address the target assignment problem 
when each agent has knowledge of all target positions, and 
a limited communication range r > 0. We introduce a class 
of distributed algorithms, called assignment-based motion, 
which provide a natural approach to the problem. Following 
the recent interest in determining the time complexity of 
distributed algorithms for robotic networks, as in [7] and 
[8], we study the worst-case asymptotic performance of the 
assignment-based motion class as the environment grows 
with n. We show that for a d-dimensional cube environment, 
[0,^(n)] d , d > 1, if the side length £(n) grows at a rate of 
at least (1 + e)rn 1 / d , where e > 0, then the completion time 
is in Q(n^ d ~ 1 ^ d £{Ti)), for all algorithms in this class. 

In Section [V] we introduce a novel control and com- 
munication algorithm, called ETSP ASSIGNMENT. In this 
algorithm, each agent computes an ETSP tour through the 
n targets, turning the cloud of target points into an ordered 
ring. Agents then move along the ring, looking for the next 
available target. When agents communicate, they exchange 
messages of O(logn) size, containing information on how 
far it is to the next available target along the ring. In 
Section IV-A1 we verify the correctness of this algorithm for 
any communication graph which contains, as a subgraph, 
the r-disk graph. In Section IV-B1 we show that when 
P( n ) > (l+e)rn 1 / d , among all algorithms in the assignment- 
based motion class, the ETSP ASSIGNMENT algorithm is 
asymptotically optimal (i.e., a constant factor approximation 
of the optimal). Finally, in Section [V-El we note that ETSP 
Assignment solves the target assignment problem even 
when there are n agents and m targets, n ^ m. 



II. Background 

In this section we introduce notation and review some 
relevant results in combinatorial optimization. 

A. Notation 

We let R denote the set of real numbers, M>o denote 
the set of positive real numbers, and N denote the set of 
positive integers. For a set finite A we let \A\ denote its 
cardinality. For two functions /, g : N — > K.>o, we write 
f(n) G 0(g) (respectively, f(n) G ^(<?)) if there exist 
JV e N and c £ K >0 such that f(n) < cg(n) for all 
n > N (respectively, f(n) > cg(n) for all n > N). If 
f(n) € 0(g) and f(n) G fl(g) we say f(n) G 6(g). 
Finally, we use the notation (mod n) to denote arithmetic 
performed modulo n G N. Thus, for an integer n G N 
we have n + 1 = 1 (mod n) and = n (mod n), and 
{n — 1, n, n + 1} = {n — 1, n, 1} (mod n). 

B. 77ie assignment problem 

Following [4], the classical assignment problem can be 
described as follows. Consider n persons who wish to divide 
themselves among n objects. For each person i, there is 
a nonempty set QM of objects that i can be assigned to, 
and cost cy > associated to each object j G QyK An 
assignment S is a set of person-object pairs (i, j) such that 
j G Q^ for all (i,j) G S. For each person i (likewise, 
object j), there is at most one pair (i,j) G S. We call the 
assignment complete if it contains n pairs. The goal is to 
find the complete assignment which minimizes Y^u j)^s Ci r 

Let Xij be a set of variables for i and j ml := {1, . . . , n}. 
For an assignment S, we write iy = 1 if (i, j) G S, and 
i.y = otherwise. Thus, the problem of determining the 
optimal assignment can be written as a linear program: 

n 

minimize £ij%ij, 

i=l jgQW 

subject to =1 V t G X, 

3GQW 

^ X.y =1 V j el, 
> 0. 

We cannot use linear inequalities to write the constraint that 
Xij's, attain only the values zero and one. However, it turns 
out, [4], that there always exists an optimal solution in which 
the Xij's satisfy our integer assumption. 

C. The Euclidean traveling salesperson problem 

Here we review some relevant results on the Euclidean 
traveling salesperson problem (ETSP). Let Q be a set of n 
points in a compact environment £ C M. d , d > 1, and let 
Q n be the set of all point sets Q C £ with \Q\ = n. Let 
ETSP(Q) denote the cost of the ETSP tour over the point 
set Q, i.e., the length of the shortest closed path through all 
points in Q. An important result, from [9], is that given a 



compact set £, there exists a finite constant a(£) such that, 
for all Q G Q n , 

ETSP(Q) < a(£)n^ d - l) ' d . (1) 

In fact, we have that in the worst-case setting, the ETSP(Q) 
belongs to Q(n { ~ d - 1 ^ d ). 

In our application of these results it will be useful to 
consider the case where the environment grows with the 
number of points. That is, we are interested in environments 
which are cubes, [0,£(n)] rf , d > 1, where £(n) is the side 
length of the cube. Applying a simple scaling argument to 
the result in (Q~|i, we arrive at the following corollary. 

Corollary 2.1 (ETSP tour length): Consider an environ- 
ment £ = [0,£(n)] d , where d > 1. For every point set 
2 G Q n , 

ETSP(Q) G 0(n'- d - 1 ^ d e(n)). 
The problem of computing an optimal tour is known to 
be NP-complete. However, there exist heuristics which can 
be computed efficiently and give a constant factor approxi- 
mation to the optimal tour. The best known approximation 
algorithm is due to Christofides [10]. The Christofides' 
algorithm computes a tour that is no more than 3/2 times 
longer than the optimal. It runs in time 0(n 3 ). Another 
method, known as the double-tree algorithm, produces tours 
that are no longer than twice the optimal, in run time 0(n 2 ). 

III. Problem formulation 

To describe the target assignment problem formally, con- 
sider ?i agents in an environment £ (n) C R d , d > 1. The 
environment £ (n) is compact for each n but may grow with 
the number of agents. For ease of presentation let £ := 
[0,£(n)] d , where £(n) > (that is, £ is a d-dimensional cube 
with side length t(n)). Each agent has a unique identifier 
(UID) taken from the set Ijjid C N. For simplicity, we 
assume that Iuid := I = {l,...,n}. However, each agent 
does not know the set of UIDs being used (i.e., agent n 
does not know it has the largest UID). Agent i G I has 
position pW G £. Two agents, i and k in T, are able to 
communicate if and only if ||pW — p' fe ]|| < r, where r > 
is called the communication range. We refer to the graph 
representing the communication links as the r-disk graph. 
Agent is kinematic model is pM = uW, where uW is a 
velocity control input bounded by v > 0. We assume that the 
agents move in continuous time and communicate according 
to a discrete time communication schedule consisting of an 
increasing sequence of time instants with no accumulation 
points, {^} fceN - We assume that \t k+1 - t k \ < t max , for 
all k G N, where t max G M>o- At each communication 
round, agents can exchange messages of length O(logn). Q 
We assume that communication round k occurs at time t k , 
and that all messages are sent and received instantaneously at 
t k - Motion then occurs from t k until t k +i- It should be noted 
that in this setup we are emphasizing the time complexity 
due to the motion of the agents. 

'The number of bits required to represent an ID, unique among n agents, 
grows with the logarithm of n. 



Let Q := {qi, . . . , q„} be a set of distinct target locations, 
qj £ £ for each j £ I. Agent i is equipped with memory 
MM, of size In this memory, agent i stores a set of 

target positions, QH C Q. These are the targets to which 
agent i can be assig ned. We let QW(0) denote agent i s 
initial target set. These positions may be known through GPS 
coordinates, or through a map of the environment. 

In this paper we assume that each agent knows the position 
of every target. That is, QM (0) = Q for each i £ X. We refer 
to this as the full knowledge assumption. To store this amount 
of information we must assume that the size of each agents' 
memory, |MW|, grows linearly with n. Our goal is to solve 
the full knowledge target assignment problem: 

Determine a control and communication law for 
n £ N agents, with attributes as described above, 
satisfying the following requirement. There exists 
a time T > such that for every agent i £ I, there 
is a unique target q 3i £ QW (0) with pM (t) = q 3i 
for all time t > T, where ji ~ jk if and only if 
i = k. 

In the remainder of the paper, we will refer to this as the 
target assignment problem. 

Remark 3.1 (Consistent knowledge): A more general as- 
sumption on the initial target sets, (0), which still ensures 
the existence of a complete assignment, is the consistent 
knowledge assumption: For each K CI, lUfcg^Q^O)] > 
/C|. In fact, it was proved by Frobenius, 1917, and Hall, 
1935 that this is the necessary and sufficient condition for 
the existence of a complete assignment [1]. • 
In the full knowledge assumption, each agent knows the 
position of all targets in Q. These positions will be stored 
in an array within each agents memory, rather than as 
an unordered set. To represent this, we replace the target 
set Q with the target n-tuple q :— (qi,...,q„), and the 
local target set QW with the n-tuple qW. Thus, in the full 
knowledge assumption, ql l l(0) := q for each i £ X. (It is 
possible that the order of the targets in the local sets qM may 
initially be different. However, given a set of distinct points 
in R d , it is always possible to create a unique ordering.) 

IV. Assignment-based algorithms with lower 

BOUND ANALYSIS 

In this section we introduce and analyze a class of deter- 
ministic algorithms for the target assignment problem. 

A. The assignment-based motion class 

The initialization, motion, and communication for each 
algorithm in the assignment-based motion class have the 
following attributes: 

Initialization: In this class of algorithms agent i initially 
selects the closest target in qW, and sets the variable currW 
(agent i's current target), to the index of that target. 

Motion: Agent i moves toward the target curr^ at 
speed v: 

^llq^-pWll' lfq cun-.^P ' (2) 

0, otherwise, 



where v > is a constant. 

Communication: As agent i communicates with other 
agents, it updates the tuple q™ "removing" targets which 
are assigned to other agents. If agent i must change curr^, 
it selects a new target in qW, that has not been removed. 
This is described more formally in the following. 



B. Lower bound on task complexity 

In order to classify the time complexity of the assignment- 
based motion class of algorithms, we introduce a few useful 
definitions. We say that agent i £ X is assigned to target 
qj* , j £ X, when currW = j. In this case, we also say target 
j is assigned to agent i. We say that agent i £ X enters 
a conflict over the target curr^, when agent i receives a 
message, msg^'l, with curr^l = curr^. Agent i loses the 
conflict if agent i is farther from curr^ than agent k, and 
wins the conflict if agent i is closer to curr™ than agent k, 
where ties are broken by comparing UIDs. 

Now we show that if agent i is assigned to the same target 
as another agent, it will enter a conflict in finite time. 

Lemma 4.1 (Conflict in finite time): Consider any com- 
munication range r > 0, and any fixed number of agents 
?i £ N. If, for two agents i and k, currW = curr^l at some 
time t\ > 0, then agent i (and likewise, agent k) will enter 
a conflict over currH in finite time. 

Proof: For each n the region £ is compact, and the 
motion for each agent is given by (O. Hence, agent i will 
reach currW in no more than di&m(£)/v time units, as will 
agent k. The condition ||pM — p' fe '|| < r will be satisfied 
within diam(£ (n))/v time units. At the next communication 
round, agent i will enter a conflict over currW. ■ 

With these definitions we give a lower bound on the 
time complexity of the task assignment problem when the 
environment grows with the number of agents. 

Theorem 4.2 (Time complexity of target assignment): 
Consider n agents, with communication range r > 0, in an 
environment £ = [0,£{n)] d , d > 1. If i(n) > (1 + e)rn^ d , 
where e £ M>o, then for all algorithms in the assignment- 
based motion class, the time complexity of the target 
assignment problem is in Sl(n( d_1 ^ d ^(n)). 

Proof: We will construct a set of target positions 
and a set of initial agent positions for which the bound 
holds. The environment £ is the d-cube, [0,£(n)] d . Divide 
the cube £ into [\n 1 / d \) d cubes, each with side length 
•£(n)/|Vi 1 / d ], and place a target at the center of each of 
the cubes until you run out. This is shown in Fig. Q] 
Notice that the distance between any two targets is lower 



Communication round for agent i. 

1: Broadcast a message, msgM, based on qW and containing curr' 1 ! 
and the UID i. 

2: Receive msgM from each agent k within communication range. 
3: for all msgM received do 

4: Based on msgM, (possibly) remove assigned targets from qW. 
5: if currM = currM then 

6: If agent i is farther from currM than agent fc, or if they are 
the same distance but i < k, remove the target given by 
currW from qW. 

7: Set curr' 1 ' to a target in qW (i.e., a target that has not been 
removed). 
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Fig. 1. Targets and agents placed on a lattice for the proof of Theorem |4.2| 
The lattice is split into 4 blocks, each containing 3 2 = 9 agents. The center 
agent of each block is shown along with its communication radius r. The 
distance between these center agents is lower bounded by ^(n)/ [n 1 / 2 ] . 
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Fig. 2. The map tour, creating an ETSP tour of seven targets. 





prev' f ' = 6 



Fig. 3. The initialization for agent i. 



bounded by £{n)/\n 1 / d ~\, and that, for sufficiently large n, 
> (1 + e)rn 1 / d /[n 1 /<*] > r - 
Next, place agent 2 at q2, agent 3 at q3 and so on so 
that pW = q;, for all i G {2, . . . , n}. From the initialization, 
we have that curr^ = i for each i G {2, . . . , n}. Now, if 
we place agent 1 in £ \ {qi, . . . , q„}, it will lose a conflict 
over any of the n — 1 occupied targets q2, . . . , q„. Thus, 
the assignment will not be complete until agent 1 reaches 
target q x . Since the distance between targets is greater than 
r, communication between agents i and k is not possible for 
any i, k G {2, . . . , n}. So, agent i G {2, . . . , n} will com- 
municate only with agent 1. Thus, prior to communication 
with agent 1, each agent i 6 {2, . . . , n} will have qM = q. 
The first time agent 1 comes within distance r of a target 
j G {2, . . . , n}, in the best-case, agent 1 will remove target 
j from qW. Now, for any deterministic method of selecting 
currH, we can place agent 1 in £ \ {qi, . . . , q„} such that 
target q x is the last target for which agent 1 will come within 
distance r. Therefore, agent 1 must come within distance r 
of each of the n — 1 assigned targets, before finally arriving 
at q x . 

Now we will lower bound the distance traveled by agent 1 . 
To do this, split the large <i-cube into [n/3 d \ smaller <i-cubes, 
or blocks, where each block contains 3 d targets. An example 
is shown in Fig. [TJ There is one target at the center of each 
of these blocks, and agent 1 must come within distance r of 
it. The distance between the center target of each block is 
lower bounded by the distance between targets, £(n) / \n 1 / d ~\ . 
Agent 1 must travel this distance at least [n/3 d \ — 1 times. 
So we have 



) ^ g n(„<^(n)). 



Path length > 

Hence, the path length is in fi(n( d-1 )/ d £(n,)). Since v G 



M>o, the time complexity is also in ^{n^ d ~ 1 ^/ d t(n)). ■ 
Remark 4.3 (£(n) < £ cr it)-' We have lower bounded the 
time complexity when £(n) grows faster than some critical 
value, £ cr it = rn 1 / d . This same type of bound appears in 
percolation theory and the study of random geometric graphs, 
where it is referred to as the thermodynamic limit [11]. 
When £(n) < £ cr it, congestion issues in both motion and 
communication become more prevalent, and a more complex 
communication and motion model would ideally be used. • 
In the next section we introduce an asymptotically optimal 
algorithm in the assignment-based motion class. 

V. The ETSP Assignment Algorithm 

In this section we introduce the ETSP ASSIGNMENT 
algorithm — an algorithm within the assignment-based mo- 
tion class. We will show that when £(n) grows more quickly 
than a critical value, this algorithm is asymptotically optimal. 
The algorithm can be described as follows. 

For each i G I, agent i computes a constant factor 
approximation of the optimal ETSP tour of the n targets in 
qM, denoted tour(qH). We can think of tour as a map which 
reorders the indices of qM; tour(qW) = (q^j-^, • ■ ■ , q^j„))i 
where a : T — > T is a bijection. Notice that this map is 
independent of i since all agents use the same method. An 
example is shown in Fig. [2] Agent i then replaces its n- 
tuple qW with tour(qW). Next, agent i computes the index 
of the closest target in qW, and calls it currW. Agent i also 
maintains the index of the next target in the tour which may 
be available, next^, and first target in the tour before currM 
which may be available, prevM. Thus, next^ is initialized 
to currW + 1 (mod n) and prevW to curr^l — 1 (mod n). 
This is depicted in Fig. [3] In order to "remove" assigned 
targets from the tuple qW, agent i also maintains the n- 
tuple, status W. Letting status ^(j) denote the jth entry in the 



n-tuple, the entries are given by 

q if agent i knows q^ J is assigned 

' to another agent, (3) 
1, otherwise. 

Thus, status W is initialized as the n-tuple (1, . . . , 1). The 
initialization is summarized in Table U 

TABLE I 

The initialization procedure for agent i. 



Initialization for agent i. 
Assumes: qH — q for each i 6 X. 

1: Compute a TSP tour of qH, tour(qH), and set qM := tour(qM ). 

2: Compute the closest target in qH , and set currM equal to its index: 
currM := argminjg;r{||q^ - pM||}. 

3: Set nextl'l := currM + 1 (mod n). 

4: Set prev' 1 ' := curr' 1 ' — 1 (mod n). 

5: Set status := X n (i.e., an n-tuple containing n ones). 



At each communication round agent i executes the algo- 
rithm COMM-RD displayed in Table IIT1 at the end of the paper. 
The following is an informal description. 



Fig.|4]gives an example of COMM-RD resolving a conflict 
between agents i and k, over curr^ = currf fe J . In this figure, 
all other agents are omitted. 

We are now ready to define the algorithm ETSP ASSIGN- 
MENT for solving the target assignment problem. 

Definition 5.1 (ETSP ASSIGNMENT): The ETSP AS- 
SIGNMENT algorithm is the triplet consisting of the initial- 
ization of each agent (see Table H), the motion law in (O, 
and COMM-RD (see Table HU, which is executed at each 
communication round. 

A. Correctness of ETSP ASSIGNMENT 

We will now prove the correctness of ETSP ASSIGN- 
MENT. It should be noted that this result is valid for any 



next'*' = next'*' = 1 




(a) Setup before the conflict over target 7. 




(b) Setup after resolution of the conflict. 

Fig. 4. The resolution of a conflict between agents i and k over target 7. 
Since agent i is closer to target 7 than agent k, agent % wins the conflict. 



communication graph which contains the r-disk graph as a 
subgraph. In order to prove correctness, let us first present 
some properties of the algorithm. 

Lemma 5.2 (ETSP ASSIGNMENT properties): During an 
execution of ETSP ASSIGNMENT the following statements 
hold: 

(i) Once target j 6 X, is assigned to some agent, the 
assignment may change, but target j remains assigned 
for all time. 

(ii) Agent i is assigned to the target curr^ which satisfies 
status^ 1 (curr [?1 ) = \ 

(iii) For agent i, statusW(j) = 0, for each j £ {prevM + 
l,prevH + 2, . . . ,nextH - 1} \ {curr^} (mod n). 

(iv) For agent i, status W (j) =0 only if target j is assigned 
to some agent k ^ i. 

(v) If, for agent i, statusW(j) = at some time t\, then 
status = for all t > ti. 

(vi) If agent i receives msgl fc l during a communication 
round, agent i will set status = for each j g 
{prev^l + 1, . . . , nextl fc l - 1} \ {curr^} (mod n). 

Proof: Statements |(ii)| and |(v)| and |(vi)| follow directly 
from the initialization and the algorithm COMM-RD. 

To see[(i)J consider an agent i, who is assigned to target j. 
Agent i's assignment can change only if it loses a conflict 
over target j. In every conflict there is a winner and the 
winner remains assigned to target j. 

Statement |(iii)| is initially satisfied since prevM + 1 = 
curr W = next'*' — 1 implies that {prevM + 1, . . . ,nextW — 
1} \ {currW} = 0. Assume that statement |(iii)| is satisfied 
before the execution of COMM-RD. At the end of COMM-RD, 
p rev W is updated to the first target before currW in the tour 



Informal description of comm-rd for agent i 

Assumes: status' 1 ' (s) = for each s 6 {prev' 1 ' + l,prevM + 
2, . . . , next ' - 1} \ {currM} ( mo d n). 

1: Broadcast msgM, consisting of the target indices, prevW, curr' 1 ', and 
nextM, the UID i, and the distance to the current target, dist' 1 '. 

2: for all messages, msgM i received do 

3: Set status' l '(j) to assigned ('0') for each target j from prevW + 

1 (mod n) to nextM — 1 (mod n) not equal to curi'M . 
4: if prevW = nextM = currM ^ curr' 1 ' then 
5: Set the status of curi'M to (because it was missed in the 
previous step). 

6: if curr' 1 ' = curr' fc ' but agent i is farther from currM than agent k 

(ties broken with UIDs) then 
7: Set the status of curi'H to assigned ('0'). 
8: if currM = curr'' 1 ' and agent i is closer than agent k then 
9: Leave currM unchanged. However, agent k will set currM to 
a new target. This target will be at least as far along the tour 
as the farther of nextM and nextM. So, set the status of nextM 
and nextM to assigned ('0'). 
10: if the status of every target is assigned ('0') then 
11: Exit ETSP ASSIGNMENT and stop motion. (This can occur only 
if there are more agents than targets and every target is assigned.) 
12: else 

13: Update curr' 1 ' to the next target in the tour with status available 
('1'), nextM to the next available target in the tour after currM, 
and prev' 1 ' to the first available target in the tour before currM. 



with status available (T). If status W (currM) = 1 then currW 
remains unchanged. If status^ (curr^) = then currM is 
increased to the first target with status available (T). Finally, 
nextM is set to the first target after currM which is available. 
Thus, at the end of COMM-RD the status of prevW, currW 
and next^ are available, and status^' (j) = for each target 
j G {prevW + 1, . . . , next^ — 1} \ {currW} (mod n). 

Statement |(iv)| is also initially satisfied since status^ = 1 ?1 
for each i el Assume Statement |(iv)| is satisfied before the 
execution of COMM-RD and that during this communication 
round agent i changes the status of a target j to assigned 
('0')- We will show that Statement |(iv)| is still satisfied 
upon completion of the execution of COMM-RD. In order 
for status to be changed, agent i must have received 
a message, msg^l, for which one of the following cases 
is satisfied: (1) Target j ^ currW lies between prev^l and 
next[ fc l on the tour; (2) There is a conflict between agents 
i and k over target j which agent i loses; or, (3) There is 
a conflict between agents i and k which agent i wins and 
next^l = j or next^ = j. 

In Case (1) either status^! (j) = or curr^l = j, and thus 
target j is assigned. In Case (2) agent k won the conflict 
implying curr^'l = j entering the communication round. 
Thus after the communication round, curr^ ^ j and target j 
is assigned to another agent. In Case (3), currW = curr^l ^ 
j, and agent k loses the conflict. In this case, agent k will 
change curr^l to the next available target on its tour. All 
targets from prev^l + 1 to next^ — 1 have been assigned. 
Also, during the communication round, agent k will receive 
msgM and determine that all targets from prevM + 1 to 
nextH - 1 are assigned. Thus, the next available target is 
at least as far along the tour as the farther of next™ and 
next^l. Thus, after the communication round, both nextM 
and next^l are assigned. ■ 

With these properties we are now ready to present the 
main result of this section. 

Theorem 5.3 (Correctness of ETSP ASSIGNMENT): For 
any fixed n G N, ETSP ASSIGNMENT solves the target 
assignment problem. 

Proof: Assume by way of contradiction that at some 
time ti > there are J G {1, . . . , n — 1} targets unassigned, 
and for all time t > t\, J targets remain unassigned. By 
Lemma I5.2l[(i)| the n — J assigned targets remain assigned 
for all time, and thus it must be the same J targets which 
remain unassigned for all t > t\. Let J denote the index 
set of the J unassigned targets. From our assumption, and 
by Lemma I5.2l|(iv)| for every t >t\ and for every i £ I, 
status W (j) = 1 for each j G J . Now, among the n — J 
assigned targets, there is at least one target to which two or 
more agents are assigned. Consider one such target, call it 
ji, and consider an agent i\ with curri 11 ) = j\. By Lemma 
14.11 agent i\ will enter a conflict over j\ in finite time. Let 
us follow the loser of this conflict. The losing agent, call 
it %2, will set status^ 2 ! (ji) = 0, and will move to the next 
target in the tour it believes may be available, call it j'2. Now, 
we know j% is not in J, for if it were J — 1 targets would 
be unassigned contradicting our assumption. Moreover, by 



Lemma l572l|(ii)| j% ^ j\. Thus, agent i 2 will enter a conflict 
over j2 in finite time. After this conflict, the losing agent, call 
it 13, will set status I* 3 ' (j'2) = (because it lost the conflict), 
and from Lemma l5^2l|(vi)| status!* 3 ! (ji) = 0. Again, agent 
i3's next target, j% must not be in J , for if it were we would 
have a contradiction. Thus, repeating this argument n — J 
times we have that agent loses a conflict over j n -j- 

After this conflict, we have status ^"~ J ^(jk) = for each k G 
{1, . . . , n — J}, where jk 1 = jk 2 if an d onr y if hi = k%. In 
other words, agent i n -j knows that all n — J assigned targets 
have indeed been assigned. Also, by our initial assumption, 
status [l ™- J l(j) = 1 for each j G J. By Lemma 15^1 [(n)] 
agent i n -/s new current target must have status available 
('1'). Therefore, it must be that agent i n -j will set cun^ n ~ J } 
to a target in J . Thus, after a finite amount of time, J — 1 
targets are unassigned, a contradiction. ■ 

The following remark displays that the ETSP ASSIGN- 
MENT algorithm does not solve the target assignment under 
the consistent knowledge assumption. 

Remark 5.4 (Consistent knowledge: cont'd): Consider as 
in Remark [3~T1 the consistent knowledge assumption for each 
agent's target set. Specifically, consider two agents, 1 and 2, 
with initial target sets Q [1 ](0) = {q 2 }, Q [2l (0) = {qi,q2}, 
and any initial positions such that p^^O) = q2, We will 
have currM = currl-'l = 2. However, agent 2 will win the 
conflict over target 2. Thus, agent 1 will set status W (2) = 0, 
and a complete assignment will not be possible. • 

B. Time complexity for ETSP ASSIGNMENT 

In this section we will give an upper bound on the time 
complexity for ETSP ASSIGNMENT. We will show that 
when £(n) > (1 + t)rn Y l d , for some e G R> , ETSP 
Assignment is asymptotically optimal among algorithms 
in the assignment-based motion class. Before doing this, let 
us first comment on the lower bound when the environment 
grows at a slower rate. 

In what follows we show that if an agent arrives and 
remains at its assigned target for sufficiently long time, then 
it stays there for all subsequent times. 

Lemma 5.5: Consider n agents executing ETSP ASSIGN- 
MENT with communication range r > and assume the 
time delay between communication rounds, t max , satisfies 
tmax < r/v. If there exists a time t\ and an agent i 
such that p^(i) = currW for all t G [h,ti + t max ], then 
pM(t) = cun-M for all t > t± + t max . 

Proof: Consider agent i, who has been at target 
currW during the entire time interval + t max ]. By 

the definition of t max there was a communication round 
at some time ti G +t max [. Agent i must have won 

any conflicts it entered during this communication round 
since we have assumed that pW(ti + t max ) = currW. Thus 
every agent k within distance r of curr^ will have set 
status^! (currW) = 0. After the communication round at £2, 
every agent k with curr^l = curr^ must be a distance greater 
than r from curr^. Since t max < r/v, any agent k that 
enters a conflict with agent i at time t > t 2 , will be at a 
distance dist^ G ]0,r[ from curr^l. Agent k will lose the 



conflict since dist^ > = dist^. Thus, agent i will remain 
at currW for all t > t\ + t rnax . ■ 
With this lemma we are now able to provide an upper bound 
on the time complexity of our scheme. 

Theorem 5.6 (Time complexity for ETSP ASSIGNMENT): 
Consider an environment £ = [Q,£(n)] d , d > 1. If 
tmax < r/v, then ETSP ASSIGNMENT solves the 
target assignment problem with time complexity in 
0{n^ d -^IH{n) + n). If, in addition, £{n) > (1 + e)™ 1 ^, 
where e G R>o, the time complexity is in Q(rv- d ~ 1 '' d t{n)), 
and ETSP ASSIGNMENT is asymptotically optimal among 
algorithms in the assignment-based motion class. 

Proof: Consider any initial agent positions, 
(0), . . . (0), and any n-tuple of target positions, 
q. In the worst-case, some agent must travel around its 
entire ETSP tour, losing a conflict at each of the first 
n — 1 targets in the tour. By Lemma 15.51 this agent 
can spend no more than t max time units at each of the 
n — 1 targets, before losing a conflict. Since each agent's 
tour is a constant factor approximation of the optimal, 
the tour length is 0(n( d ~ 1 ^ d £(n)) (see Theorem 12. It . 
The agent will not follow the ETSP tour exactly because 
it will enter a conflict over each of the n — 1 targets 
before actually reaching the target. However, the resulting 
path is no longer than the ETSP tour (since the agent 
could just follow the ETSP tour exactly if that happened 
to be the shortest path). Hence, the time complexity is 
0( n (d-i)/d£( n ) + t max (n - 1)) G 0{n( d - 1 V d l{n) + n). 
If i{n) = (2 + e)r?i 1 / d , with e G R>o, we can combine 
this with Theorem 14.21 to get a time complexity in 

e( n (<»-i)/<^( n )). ■ 

Notice that when £(n) satisfies the bound in Theorem 15.61 
and £(n) G 0(n 1 / d ), the time complexity is in 0(n). 

We have given complexity bounds for the case when r and 
v are fixed constants, and £(n) grows with n. We allow the 
environment £(n) to grow with n so that, as more agents are 
involved in the task, their workspace is larger. An equivalent 
setup would be to consider I to be fixed, and allow r and v 
to vary inversely with the n. That is, we can introduce a set 
of parameters, £ = 1, and r(n) and v(n) such that the time 
complexity will be the same as for the parameters r, v, £(n). 

Corollary 5.7 (Scaling radius and speed): Consider 
n agents in the environment £ = [0, with speed 
v(n) := v/£(n), and communication radius f(n) := r/£(n), 
where £{n) > (1 + e)rn 1 / d , and e G R> . Then ETSP 
ASSIGNMENT solves the target assignment problem with 
time complexity in <d(n( d ~ 1 ^ d £(n)). 

Scaling the communication radius r inversely with the 
number of agents arises in the study of wireless networks 
[12]. In wireless applications there are interference and 
media access problems between agents in the network. Since 
the agents are in a compact environment, the only way to 
limit this interference is to scale the communication radius 
inversely with the number of agents. Scaling the agent speed 
inversely with n appears in the study of the vehicle routing 
problem in [7]. The inverse scaling is required to avoid 
collisions in the presence of traffic congestion. 



C. Communication and computation complexity 

In our notion of time complexity we have emphasized 
the complexity due to the motion of the agents. Here we 
will briefly classify the complexity of computation and 
communication for ETSP ASSIGNMENT, (i) Initialization: 
As reviewed in Section III-CI we can compute a constant 
factor approximation ETSP tour in time 0(n 2 ). This is 
the most expensive computation and thus the complexity 
of initialization is in 0(n 2 ). (ii) Communication complexity 
per round: At each round agent i broadcast a message of 
length O(logn), msgW, and we consider this to be one unit 
of communication. In the worst-case, each agent receives n 
messages, and so, the worst-case communication complexity 
is in 0(n) [8]. (iii) Computation complexity per round: For 
each message received, agent i sets status^ (,s) = for s 
from prevM + l to next[ fe ] — 1. In the worst-case, this operation 
is 0(n) and must be performed for n messages. This is the 
dominant computation in COMM-RD and thus the worst-case 
computation complexity in each round is 0(n 2 ). 

It should be noted that in the case when the communica- 
tion graph is not even connected (let alone complete as is 
required to achieve these worst-case bounds), the complexity 
will be considerably lower. 

D. Simulations 

We have simulated ETSP ASSIGNMENT in R 2 and R 3 . 
To compute the ETSP tour we have used the concorde 
TSP solver@ A representative simulation for 15 agents in 
[0, 100] 3 C K 3 with r = 15 and v = 1 is shown in Fig. 
[5] The initial configuration shown in Fig. |5(a)| consists of 
uniformly randomly generated target and agent positions. 

E. The case of n agents and m targets 

It should be noted that the ETSP ASSIGNMENT algorithm 
works without any modification when there are n agents and 
m targets. If m > n, at completion, n targets are assigned 
and m — n targets are not. When, m < n, at completion, 
all m targets are assigned, and the n — m unassigned agents 
come to a stop after losing a conflict at each of the m targets. 
The complexity bounds are changed as follows. 

The lower bound on the assignment-based motion class 
in Theorem 14.21 holds when m > n, and £(n) > (1 + 
e)rm 1 / d (notice the m instead of n). The bound becomes 
fi(^(n)m- 1 / d ?i). If m = Cn where C E K>i, (i.e., m > n 
but they grow at the same rate), then the bound becomes 
fi(^(n)n( d - 1 )/ d ). 

The upper bound on ETSP ASSIGNMENT holds for any 
n and m, and becomes O^njN^ 1 ^' ' d ), where N := 
min{n, m}. So our final result would be that if m = Cn 
where C G R>i and when £(n) > (1 + e)rm 1 / d , then 
ETSP Assignment solves the target assignment prob- 
lem in 9(£(n)n (d - 1)/d ). That is, among all algorithms in 
the assignment-based motion class, ETSP ASSIGNMENT is 
asymptotically optimal. 

2 The concorde TSP solver is available for research use at 

http : / / www . t sp . gatech . edu/ concorde/ index . html 




(a) Initial agent and target positions. (b) Positions after 30 time units. 




(c) Positions after 90 time units. (d) Complete target assignment. 



Fig. 5. Simulation for 15 agents, with v = 1 and r = 15. in an environment 
[0, 100] 3 . The targets are spheres and the agents are cubes. An edge is drawn 
between two agents when they are within communication range. 



VI. Conclusions 

We have developed the ETSP ASSIGNMENT algorithm for 
solving the full knowledge target assignment problem. We 
derived worst-case asymptotic bounds on the time complex- 
ity, and we showed that among a certain class of algorithms, 
ETSP Assignment is asymptotically optimal. There are 
many possible extensions of this work. We have not given a 
lower bound on the time complexity of ETSP ASSIGNMENT 
when l(n) < £ cr it- Also, the problem is unsolved under the 
more general consistent knowledge assumption. We would 
like to extend the ETSP Assignment algorithm to agents 
with nonholonomic motion constraints. Also, it would be 
interesting to consider a sensor based version of this problem, 
where agents acquire target positions through local sensing. 
Finally, to derive asymptotic time bounds, we made some 
assumptions on the communication structure at each commu- 
nication round. An interesting avenue for future study would 
be to more accurately address the communication issues in 
robotic networks. 
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TABLE II 

Communication Round (comm-rd) for agent i. 



Name: comm-rd 

Goal: Obtain information on assigned targets. 

Assumes: (i) Knowledge of the n-tuple q, and a method for 
computing a constant factor TSP tour of the n targets, 
tour, (ii) A communication range r > 0. 


1 


Compute aist L J := p L 1 — q r - q . 

ii cum 1 J 




2 


r i r 'i r 'i r i r *' 
Broadcast msg 1 ' 1 := (prev 111 , curri 11 , next 111 , i, dist w 


Km,, 


3 


Receive msg' fc l, from each k ^ i satisfying Hp' 1 ' — 


P l 'II < r. 


4 


for all msgM received do 




5 


for s = prev' fe l + 1 to next' fe ' — 1 (mod n) do 




6 


if s 7^ curr' 1 ' then 




7 


Set status := 




8 


r* i rZ"1 rzi f'l 
if prev 1 > = next 1 ' = curr 1 ' jt curr 111 then 




9 


Set status W(currM) := 




10 


if curr' 1 ' = curr' fc l then 




11 


if (distM > distW) OR (distW = distal AND 


i < k) then 


12 


Set status M (currM) := 0. 




13 


else 




14 


if nextM ^ currM then 




15 


Set statusW(nextW) := 0. 




16 


if next' fe ' ^ curr' 1 ' then 




17 


Set statusW(nextM) := 0. 




18 


if status M (j) = for every target j then 




19 


Exit ETSP ASSIGNMENT and stop motion. 




20 


else 




21 


while status'*! (curr' l l)=0 do 




22 


currM := currM + 1 (mod n). 




23 


Set next' 1 ' := curi'M + 1 (mod n). 




24 


while status'*! (next' l ')=0 do 




25 


next' 1 ! := next' 1 ' + 1 (mod n). 




26 


while status W (prevM)=0 do 




27 


prev' 1 ' := prev' 1 ' — 1 (mod n). 
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